Dynamics of Passive-Scalar Turbulence 
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We present the first study of the dynamic scaling or multiscaling of passive-scalar turbulence. 
For the Kraichnan version of passive-scalar turbulence we show analytically, in both Eulerian and 
quasi-Lagrangian frameworks, that simple dynamic scaling is obtained but with different dynamic 
exponents. By developing the multifractal model we show that dynamic multiscaling occurs in 
passive-scalar turbulence only if the advecting velocity field is itself multifractal. We substantiate 
our results by detailed numerical simulations in shell models of passive-scalar advection. 

PACS numbers: 47.27.i, 47.53.+n 



Important advances have been made over the past 
decade in understanding the statistical properties of the 
turbulence of passive scalars, e.g., pollutants, and passive 
vectors, e.g., weak magnetic fields, advected by a fluid 
If the advecting velocity is stochastic and of the Kraich- 
nan type 0, then, in some limits, it can be established 
analytically that passive-scalar and passive- vector turbu- 
lence show anomalous scaling or multiscaling of structure 
functions. These are the only turbulence problems for 
which such multiscaling can be proven analytically, so 
it is important to use them as testing grounds for new 
ideas about multiscaling in turbulence. We examine the 
dynamic-scaling properties of passive-scalar turbulence 
in the light of the recent systematization of dynamic mul- 
tiscaling in fluid turbulence Q. 

Our principal results, obtained analytically and numer- 
ically, illustrate important principles that appear, at first 
sight, to be surprising. We find, e.g., that the dynamic 
exponents depend via bridge relations only on the equal- 
time scaling exponents of the velocity field. Thus, even 
though equal-time structure functions for the passive- 
scalar problem display multiscaling, we show analytically 
that they exhibit simple dynamic scaling if the advect- 
ing velocity is of the Kraichnan type. Our study of the 
Kraichnan model yields the only analytical results ob- 
tained so far for time-dependent structure functions in 
any type of turbulence. Dynamic multiscaling is obtained 
only if the advecting velocity field is itself intermittent as 
we demonstrate numerically. 

The quest for a statistical characterization of turbu- 
lence begins most often with the equal-time, order-p, 
velocity- m structure functions Sp(£), i.e., the order-p mo- 
ments of the probability distribution functions (PDFs) 
of velocity differences at the length scale I. The equal- 
time exponents Q are defined by {I) ~ ftp , valid 
for the inertial range jjj < ^ < L, where rjd is the 
dissipation scale and L the length at which energy is 
pumped in. Kolmogorov's simple scaling 0(K41) yields 
Cp' Kil — p/3, but subsequent work [j| suggests signif- 
icant corrections for p > 3 and multiscaling with a 
nonlinear, convex, monotonically increasing function of 
p. The generalization of such multiscaling to dynamic 



multiscaling is subtle and has been elucidated only re- 
cently 0, S S 13 : It is expected that Eulerian- velocity 
time-dependent structure functions lead to trivial dy- 
namic scaling with all Eulerian (E ) dynamic exponents 
Zp- £ = 1. Nontrivial dynamic exponents can be ob- 
tained 01 from dynamic-multiscaling ansatze of the form 
Tp ~ t z p , where the times are extracted from La- 
grangian or quasi-Lagrangian time-dependent structure 
functions. A generalization of the multifractal formal- 
ism |5| to the case of time-dependent structure func- 
tions U 0, Hi shows that dynamic and equal-time mul- 
tiscaling exponents must be related by bridge relations; 
and these bridge relations depend crucially Q on how 
time scales are extracted: In particular, time scales ex- 
tracted via time derivatives of time-dependent structure 
functions are different from those obtained from integrals. 

The advection-diffusion equation for the passive scalar 
field 9(x,t) at point x and time t is 



d t + UidiO = K,d it + f s , 



(1) 



where k is the passive-scalar diffusivity and fg an exter- 
nal force. The advecting velocity u should be obtained 
from solutions of the Navier-Stokes equation, but, to in- 
vestigate equal-time multiscaling of passive-scalar struc- 
ture functions, it has proved fruitful to use the Kraichnan 
ensemble in which each component of u is a zero-mean, 
delta-correlated Gaussian random variable with 



( Ui (x,t) Uj {x + r,t')) = 2D lJ {v)8{t-t'). 
The Fourier transform of has the form 



Dij(q) cx q 
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(2) 



(3) 



where q is the wave- vector, d the spatial dimension, 77 the 
dissipation scale, L the large forcing scale, and £ a param- 
eter. The term inside square brackets assures incompress- 
ibility. In real space, if we write Dij(r) = D°Sij — |dy (r) 
and take the limits L — > 00 and n 



£>ir C 



0, we get 



(4) 
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with Di a normalization constant. For < £ < 2 
this model shows multiscaling of equal-time, order-p, 
passive-scalar structure functions . The constant D° = 
2 Jp 00 Dij(q)d d q oc O(L^) diverges in the limit L — ► oo. 
The external force /g is also a zero-mean, Gaussian ran- 
dom variable which is white-in-time with the variance 
(/e(x,t)/ fl (y s t')> = C (^) *(t - O. where the l™> 
tion C(x/L) is confined to large length scales. Moreover, 
/# and u are statistically independent. 

The quasi-Lagrangian transformation of any Eulerian 
field tp(x,t) is defined by $(x, i) = ^[x + R(i; tq, 0), t], 
where R(i;r ,0) is the Lagrangian trajectory passing 
through r at t = 01- Thus the quasi-Lagrangian 
version of Eq. (|TJ is 



ft. 



•fii(0)] di§ = KdJ + fg. 



(5) 



The equal-time multiscaling properties of passive-scalar 
structure functions in the Kraichnan model, well under- 
stood in the Eulerian framework, remain unchanged in 
the quasi-Lagrangian framework [l| . We concentrate on 
the time-dependent structure functions 



J*(r, fa,...,^}) = <[<ty(x,ii,r).. 



M P ,r)]), (6) 



where is or in the Eulerian or quasi-Lagrangian 
frameworks, respectively, the angular brackets denote an 
average over the PDFs of u and fg, and 50(x,£,r) = 
0(x + r, t) — </>(x, i). Dimensional analysis yields the in- 
dependent characteristic time T{r) ~ gfr^ ~ r 2 ~^, 

whence = 2 — £ for all p. Like the K41 result 
Zp = 2/3 for fluid turbulence, we expect this pre- 
diction to be valid for the passive-scalar case in La- 
grangian or quasi-Lagrangian frameworks. For p = 
2, jf(r,t) = 2C+(0,t) -2C*(r,f), where C+(r,t) = 
(cf)(x + r, i)0(x, 0)). Equations for C*(r, i), obtained 
from Eqs. (JTJ and © by Gaussian averaging over 
the PDFs of u and fg are (in the n — » limit of 

relevance to turbulence) 

ftC e (r,<) = D°(L)d u C e ~ £%C e ; (7) 
ftC^(r,i) = (l>% - D y ) ft^C* ~ dij^dijC 6 . (8) 

Spatial Fourier transforms of Eqs. J7J or J5J) then yield 
C^(ci,t) ~ exp[— t/r^(q)], with characteristic time scales 
r e = [Z)°(X)<j]~ 2 [Eq. (0] and t<? ~ q^ 2 [Eq. ©]. The 
last term of Eq. J7J) diverges as L — > oo, a signature of the 
sweeping effect; however, for a fixed integral scale L, the 
Fourier transform of Eq. Q implies z\ = 2 (cf. = 1). 
This sweeping divergence is removed in Eq. JHJl, whose 
Fourier transform, in the limit L — > oo, gives zf = 2 — £, 
which agrees with the dimensional prediction. These re- 
sults are remarkable for two reasons: (1) this is the first 
calculation of dynamic scaling exponents in any form of 
turbulence and (2) this shows explicitly how sweeping 



effects are removed in the quasi-Lagrangian representa- 
tion. Similar, but more cumbersome, calculations yield 
zf = 2 and z% = 2 — £, which we substantiate by our 
numerical studies below [l5j ]. 

A turbulent velocity field obeying the Navier-Stokes 
equation does not have simple statisticalproperties as in 
Eq. (0) , so we use the multifractal model |{| and its exten- 
sion to dynamic multiscaling [j| . To eliminate the sweep- 
ing effect, we consider the structure functions, namely 
Fp(r,t), namely, 



K I d^iLf+^-^gv^-L), (9 ) 
[9 l )p J i L t pJi 



where 9 is assumed to possess a range of universal scal- 
ing exponents h G X = (h rn i n , h max ). For each h in 
this range, there exists a set C R 3 of fractal di- 
mension D(h), such that S9 ^ r ^ cx {j^) h for x 6 S/j , 

with 9l the passive-scalar variable at the forcing scale 
L. G p,h {^r—) has a characteristic decay time T p> h ~ 
r/5u(r) £P' h (0) = 1, and the velocity field is 

also multifractal with a range of universal scaling expo- 
nents g £ 1 = {g mm ,9max)- Following Ref. 3f we de- 
fine the order-p, degree-M, integral time scale T p '^(r) = 

(l/M) 

f °°F»(r,t)tW-Vdt 



scale T 



p,M 



and the derivative time 



gpj-JT p (r, i) 



t=o 



with cor- 



responding dynamic-multiscaling exponents defined via 



9,D 



r z p,M and T p ° M 



s,r> 

y Z p,M 



To calculate 



d^{h){jf+vh-D(h) 

% L Jo 



Q p ' h (—)dt, 



Tph 



(10) 

e.g., we substitute the multifractal form ©, do the 
time integral first, and use the scaling ansatz for t Pj h 



and we obtain T p y (r) oc r 1_c * (((^((fa)" 1 ). Following 
Ref. ^| E3l ! if we now assume that the dominant contri- 
bution to 



(W(*u)-*) » (W>((5«)-«) 



(11) 



we have z. 



p,M 



1 — IC-il- Similar calculations yield the 



following more general bridge-relations 



p,m 



M ' 



9,1 
Z p,M 



le 



Ml 



(12) 



which do not depend on p. However, this does not mean 
we have simple dynamic scaling. For a velocity field 
which multiscales, i.e., for which C, M /M ^ £™, we have 
dynamic multiscaling. For the Kraichnan model a similar 



multifractal formalism predicts zl'^ = 2 — £ 
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We now substantiate our analytical and multifrac- 
tal results by detailed numerical simulations of shell- 
model analogs of the two passive-scalar problems dis- 
cussed above. For this purpose we consider two very 
similar shell models, A and B, of the general form |lllll2| 

[j t + nk m ]e m {t) = /<i>2,;i + <W/(*)- (13) 

For model A, which corresponds to the Kraich- 
nan model, [( W 2 ) (^m+i^m-i ~ 

C-l^m+l) + {-km-l/2){6* m _ x U* m _ 2 + e* m _ 2 U m ^x) + 

(k m+1 /2)(6* m+2 u m+ i + 0* n+1 u* n+2 )}. Here the asterisk 
denotes complex conjugation, 9 m and u m are, respec- 
tively, shell-model analogs of the Fourier components of 
the passive scalar and velocity in the quasi-Lagrangian 
framework, k m — 2 m ko and ko = 1/16. The shell 
velocity is a zero-mean, Gaussian random complex vari- 
able with covariance (w m (t)u* (t')) = D m S min 5(t — £'), 
D m = , f(t) is random, Gaussian and white-in- 
time and independent of u m . In model B, $> B - = 

ra.Ou 

[k m (0m+lU m ~l — ^m-l^m+l) — (^m-l/2) {@m-lU m - 2 + 

i)(-k m+ i/2)(6 m+2 u m+1 + m+1 u m+2 )]* . This 
is the shell-model analog of a passive scalar advected by 
a Navier-Stokes velocity field. Here the shell velocity 
obeys the GOY shell-model equation 0, ^| 

+ vk 2 m ]u m = iV m ^u + <W/™, (14) 

where T„ hilil = [k m u m+ iu m+2 - 5fc m _ 1 u m _iu m+1 - (1 - 
5)k m — 2 u m -xU m - 2 \* with 5 = 1/2. For both models 
A and B, u_i = uq = 6-\ = 6$ = Q, and the cou- 
plings are limited to next-nearest-neighbor shells, hence 
the dynamic scaling or multiscaling properties of time- 
dependent structure functions should be akin to those 
of quasi-Lagrangian variables Q (so we use 9 m and 

U m ). Furthermore, when v — > and f u — > 0, E e = 

J2 m =i I ® m I 2 i s conserved, where TV is the total number 
of shells. We use a weak, order-one, Euler scheme asso- 
ciated with the Ito formulation of Eq. 1|13|) to integrate 
model A and a second-order Adams-Bashforth scheme to 
integrate model B. The different parameters used in our 
simulations, e.g., the large-eddy-turnover-time T£, are 
given in Table (JU). The shell- model equal-time passive- 
scalar structure functions are S p (m) = ([9 m m ] p / 2 ^ ~ 
-C" 

km p . Model A exhibits equal-time multiscaling for 
< £ < 2; we use £ = 0.6 here. The equal-time mul- 
tiscaling exponents for model B are given in the second 
column of Table Hll Our results for the equal-time multi- 
scaling exponents ( p for both models A and B agree well 
with previous studies (Ref. f° r model A and Ref. 0] 
for model B). We now define the order-p, time-dependent, 



Model 


K 


St 




Ttr 




A 


2 -14 


2 -24 


~ 2 2i St 


5 x 10 4 Ti 


10 5 r L 


B 


5 x 1(T 7 


1Q -4 


10 s 8t 


5 x 10 4 t l 


10 5 T L 



TABLE I: The diffusivity n, the time-step St, and the box- 
size eddy turnover time tl = X/koUrms that we use in our 
simulations of models A and B. Data from the first T tr time 
steps are discarded so that transients can die down. We then 
average our data for time-dependent structure functions for 
an averaging time T av . For model A we use £ = 0.6. The 
number of shells N = 22 for both the models. 



order (p) 


& 


8,1 
z p,l 


6,D 
Z V,2 


1 


0.34 ±0.001 


0.52 ±0.03 


0.63 ±0.03 


2 


0.63 ±0.001 


0.53 ±0.03 


0.64 ±0.03 


3 


0.87 ±0.001 


0.56 ± 0.005 


0.64 ± 0.005 


4 


1.07 ±0.001 


0.56 ± 0.005 


0.64 ±0.01 


5 


1.24 ±0.004 


0.56 ± 0.005 


0.64 ±0.01 


6 


1.38 ±0.006 


0.57 ±0.007 


0.64 ±0.02 



TABLE II: Order— p (Column 1) multiscaling exponents for 
1 < p < 6 from our simulations of model B: equal-time ex- 
ponents Cp (Column 2), integral-scale dynamic-multiscaling 

exponent 2*'{ of degree-1 (Column 3) and derivative-time 

dynamic-multiscaling exponents (Column 4) from our 

calculations. We obtain (1) z e ^\ = 0.56 ± 0.005 from the 

bridge-relation fl2J and C-i = -0.44 ± 0.005; (2) z°f = 
0.645 ± 0.0001 from the bridge relation and £2 = 0.709 ± 
0.0001. See text for error estimates. 

passive-scalar structure functions for our shell models: 

F p (m,t) = — ^ ([8 m (0)d m (t)] p/2 ) ■ (15) 
For model A an analytical calculation yields [cf. Eq. J7J] 

F 2 (m,t) = S 2 (m)exp[-jfe^A(Ot], (16) 

where = [2< 2 «- 2 ) ± 2-( 2 «~ 2 )) ± (2« + 2~«) ± (2«~ 2 ) + 
2 _ ^ _2 ^], whence z 2 = 2 — £. Similar relations can be 
derived for p > 4 but the complexity increases with p. 

For model A we fit an exponential to each of F p (m, t) 
up to a time i^, such that -g7~p = /i, to extract a 
characteristic decay rate T p (m). We use fi — 0.7; we 
find that values of fi from 0.5 to 0.9 do not change 
our results significantly. Representative plots for p = 4 
are shown in Fig. a) . The slopes of these and sim- 
ilar plots yield zf = 1.400 ± 0.005, z{ = 1.400 ± 005, 
z\ = 1.40 ± 0.01, which clearly support our analytical 
result z\ = 2 - £ with £ = 0.6. For model B the 
imaginary parts of F p {m 1 t) are negligible compared to 
their real parts. Henceforth we consider only the real 
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FIG. 1: (a) Plot of the order-4 time-dependent structure func- 
tion for model A, with £ = 0.6, versus t/rL, for shells m — 6 
to 13 (the last few are not clearly visible). Exponential func- 
tions [like Eq. 116H 1 are good approximations for these time- 
dependent structure functions. The inset shows a log-log plot 
of rf (m) against the wavevector fc m ; the slope of the straight 
line (least-squares fit) gives zf ~ 1.40 ~ 2 — £. (b) A plot 
as in (a) for model B, for p = 6. The inset shows a log-log 
plot of T®2 versus the wavevector fc m ; the straight line is the 
least-squares fit to the points. 

part. We use two different sampling rates, 50 x St for 
4 < m < 8 and 10 x 5t for 9 < m < 13, respectively. 
For extracting T^ 2 we extend F p (m,t) to negative t via 
F p (m,~t) = F p (m,t) and use a centered, sixth-order, 
finite-difference scheme to find ■^■F p (m,t)\ t _ Q . A log- 
log plot of Tp 2 versus k m now yields the exponents z p ® 
given in Table (|lT|) . A comparison of our results with the 
multifractal prediction of Eg. 1121) is shown in Tabic 01. 
We have computed the integral time scale only for M = 1 , 

t — z^' 7 

i.e., T^ x {m,t u ) = J Q U F p (m,t)dt ~ k m pA 18]. In princi- 
ple we should use t u — > oo but, since it is not possible to 
obtain F p {m, t) accurately for large t, we select an upper 
cut-off t u such that F p (m,t u ) — a, where, for all m and 
p, we choose a = 0.7. We have checked that our results 
do not change if we use 0.3 < a < 0.8. The slope of a 

log-log plot of Tp tl (m) versus k m now yields shown in 
Table JD}. The difference between the dynamic-scaling 
exponents of integral and derivative type (Columns 3 and 
4 of Table JllJ), respectively) is a clear signature of dy- 
namic multiscaling. 

To obtain the error-bars on the equal-time and dy- 
namic exponents, we carry out 50 runs, each averaged 
over a time T av given in Table ©. We thus obt ain 50 dif- 
ferent values for each of these exponents. The mean val- 



ues of these 50 exponents are quoted here [Table 10]; and 
the root-mean-square deviation about the mean value 
yields the error estimates shown in Table (0 . Note that 
the low-order dynamic structure functions decay slowly 
and hence, for a fixed averaging time, low-order dynamic 
multiscaling exponents have larger error bars than those 
of slightly higher-order exponents. 

We have shown that, for passive-scalar turbulence, 
simple dynamic scaling is obtained if the velocity field is 
of the Kraichnan type. Dynamic multiscaling is obtained 
if the advecting velocity is itself multifractal. The study 
of dynamic multiscaling in passive-vector turbulence is 
more complicated because of the dynamo effect and will 
be dealt with elsewhere. The experimental study of such 
dynamic multiscaling remains an important challenge. 
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